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Abstract 

We prove the following sharp upper bound for the gradient of the Neumann 
semigroup Pt on a (i-dimensional compact domain O with boundary either C 2 - 
smooth or convex: 

HVPtHi-Kx, < t{d+1)/2 , t > 0, 

where c > is a constant depending on the domain and || • ||i_>oo is the operator 
norm from L 1 (0) to L°°(Q). This estimate implies a Gaussian type point- wise 
upper bound for the gradient of the Neumann heat kernel, which is applied to the 
study of the Hardy spaces, Riesz transforms, and regularity of solutions to the 
inhomogeneous Neumann problem on compact convex domains. 
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1 Introduction 

Let Q be a d-dimensional compact Riemannian manifold with boundary dQ. Let N be the 
inward unit normal vector field of dfl. If dQ is C 2 -smooth, then the second fundamental 
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form of dQ is defined as 



l(vi,v 2 ) := -(V Vl N,v 2 ), v ly v 2 e TdQ, 

where TdQ is the tangent space of dQ. We call the boundary dQ (or the manifold Q) 
convex if I(v, v) > for all v G TdQ. 

Let Pt be the Neumann semigroup generated by A, the Laplacian on Q with the 
Neumann boundary condition. Let p t (x,y) be the Neumann heat kernel, which is the 
density of P t w.r.t. the Riemannian volume measure. For any p, q > 1, let || • \\ p ^ q denote 
the operator norm from L P (Q) to L q (Q). When dQ is C 2 -smooth, it is easy to prove the 
following uniform gradient estimate of P t , which is important in potential analysis of the 
Neumann Laplacian as shown in Section 4 below. 

Theorem 1.1. Let Q be a compact Riemannian manifold with C 2 -smooth boundary. Then 
there exists a constant C > such that 



(1-1) llVPtlli-oo <jS^, *>0. 

Consequently, letting p be the Riemannian distance on Q, one has 



C 



p(x,y) 



2-, 



(i.2) \vp t (;y)(x)\ < ^+m exp [~ ct r t>0 ' x '^ Gfi 

for some constants C, c > 0. 

We remark that (11 .11) is sharp (for short time) since the equality holds for the classical 
heat semigroup on M. d and some constant C > 0. Since the boundary is C 2 -smooth such 
that the second fundamental form is bounded below, the above theorem can be proved 
by using the reflecting Brownian motion and exponential moments of its local time (see 
Section 2). 

When dQ is merely Lipschitzian such that the second fundamental form is not well 
defined, the argument presented in Section 2 is no longer valid. Indeed, for general 
manifolds with Lipschitzian boundary, even the existence and uniqueness of the reflecting 
Brownian motion is unknown. Nevertheless, for compact convex domains on M d , the 
reflecting Brownian motion has been constructed by Bass and Hsu in [5], [6] , which enables 
us to derive the above gradient estimates. 

Theorem 1.2. Let Q be a compact convex domain in M. d . Then and hold for 
some constants C, c > and p(x, y) — \x — y\. 

The proofs of these two theorems will be given in the next two sections respectively. In 
Section 4, we introduce some applications of our results to the study of the Hardy spaces, 
Riesz transforms associated to the Neumann Laplacian, and regularity of solutions to the 
inhomogeneous Neumann problem on compact convex domains. 
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2 Proof of Theorem 



1.1 



We first observe that for the proof of ( 11.11) it suffices to consider t £ (0, 1]. Let dx denote 
the Riemannian volume measure on Q. Since VPtf = VPtf, where / := / — J Q f(x)dx, 
we may assume that / itself has zero integral on Q. Let Ai(> 0) be the first Neumann 
eigenvalue on Q. We have 



(2.1) ||PJ 2 || 2 <e-^||/|| 2 , t>0. 

On the other hand, by taking e.g. s = t in [T9l Corollary 1.4], we obtain 

\B{x,Vt)\-\P t f{x)\<C x [ P2t\f\(y)dy<C 2 \\f\\ u t £ (0, 1] 

JB{x,<Ji) 

for some constant C\ > 0, where B(x,\/t) := {y £ f2 : p(x,y) < \/t} with volume 
\B(x, Since Q is compact, we have \B(x, y/t)\ > ct d l 2 for some constant c > and 

allt > 0, x £ Vt. Therefore, 



(2-2) ||PtU<^ «G(0,1] 

for some constant C 2 > 0. Now, if flH|) holds for t < 1, then §TQ) and Q yield 

HVPt/IU = HVPi/aP^iPi/a/IU < C 3 ||P-iPi /2 /||2 

< Cae-^-^H^/lb < C^-^H/IU, f >1 

for some constants C3, C4 > 0. Therefore, fll.ip also holds for t > 1. 

From now on, we assume that t < 1. Let K, a be two constants such that Ric > —K 
holds on Q and I> —a holds on <9f2. By [181 Theorem 1.1(7)], we have 



(2-3) 



IWI 2 ^ 71 W^i^ 2 ~ ( P */) 2 ) E / e^'-^'da, / £ C(M), 



where (Z s ) s >o is the local time on dQ for the reflecting Brownian motion on Q. Next, 
according to [T6J Proof of Lemma 2.1], there exists a constant c > such that 

Ee ala < e c{s+v ~ s \ s<l. 
Therefore, it follows from (12. 3p that 

||VP/||L<^||P/ 2 ||oo, feC(M) 
holds for some constant C5 > 0. Replacing / by P t f and using (12.21) . we arrive at 
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||VP 2t /|U< i ^||/|| 1 , feC(M) 

for some constant Cq > 0. This proves fll.ip for t < 2. 
Finally, since fll.ip is equivalent to 

C 

sup |Vpt(-,j/)(a;)| < t(d+1)/2 , t > 0, 
the inequality (11. 2p follows from the self- improvement property as in jTOj Theorem 4.9]. 



3 Proof of Theorem EH 

We shall make use of the reflecting Brownian motion determined as solutions to the 
Skorokhod equation 



(3.1) X t = x + W t + [ N(X s )dl s , t>0, 

Jo 

where x G Q, Wt is a d- dimensional Brownian motion on a complete filtrated probability 
space ($, JP t ,F), X t is a continuous adapted process on fl, and l t is a predictable contin- 
uous increasing process with l = which increases only when X t e <9f2. If (X t , Z t ) solves 
( 13 . 1 p for some (i-dimensional Brownian motion Wt, we call X t the reflecting Brownian 
motion on ^starting from x, and call It its local time on dQ. 

Lemma 3.1. For any d- dimensional Brownian motion W t and any x G Q, (\3.1\ has a 
unique solution. 

Proof, (a) Uniqueness. Let (X t ,l t ) and (X t ,l t ) be two solutions to ( 13. ip . By the Ito 
formula, 

\X t -X t \ 2 = 2 [ (X s -X s ,N(X s ))dl s + 2 [ (X s -X s ,N{X s ))dl s . 
Jo Jo 

By the convexity of dQ we see that (y — z, N(y)) < if y G dQ and z G SI. Moreover, 
since dl s = (correspondingly, dl s = 0) for X s dQ (correspondingly, X s ^ dQ), we 
conclude that \X t — X t \ 2 = for alH > 0. 

(b) Existence. By using the regularity of the associated Dirichlet form, [61 Theorem 4.4] 
ensures the existence of a reflecting Brownian motion. According to Theorem 1], this 
reflecting Brownian motion solves the Skorokhod equation ( 13. ip for some <i-dimensional 
Brownian motion Wt, i.e. there exists (Xt, k) such that X t is a continuous adapted process 
on Q, l t is a predictable continuous increasing process with Iq = which increases only 
when Xt G dQ, and 
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(3.2) x t = x + W t + [ N(X s )dl s , t>0 

Jo 

holds. We aim to prove that for any d- dimensional Brownian motion W t , the equation 
has a solution. Due to (a), this follows from the Yamada-Watanabe criterion [2U] (cf. the 
proof of [U Theorem 5.9]). For readers' convenience, we present below a brief proof. By 
the uniqueness of solutions to (13. 2p . (X,l) is determined by W; that is, there exists a 
measurable function 

F = (F u F 2 ) : C([0, oo); M d ) -> C([0, oo); Q) x C 7 ([0, oo); [0, oo)) 

for a- fields induced by the topology of locally uniform convergence, where C/([0, oo); [0, oo)) 
is the space of all non-negative continuous increasing functions on [0, oo) with initial data 
0, such that 

(X,!) = F(W) = (F 1 (W),F 2 (W)). 

So, letting fi be the Wiener measure on C([0, oo); M d ) (i.e. the distribution of W), we 
have 

F l {u) t = x + uj t + [ N{F l {u) s )dF 2 {u) s , t>0 
Jo 

for /i-a.e. uj G C([0, oo); M d ). Therefore, for any c?-dimensional Brownian motion W, which 
has the same distribution fi, 

F 1 (W) t = x + W t + [ N(F 1 (W) s )dF 2 (W) s , t>0 
Jo 

holds a.s. This means that (X, 1) := (Fi(W), F 2 (W)) solves the equation (EI]). □ 

Lemma 3.2. Let Xf be the unique solution to (\3.1\) for Xq = x G Q. For f G ^(fi), t/ie 
c/ass o/ all bounded measurable functions on Q, let Ptf(x) := u(t, x) solve the Neumann 
heat equation 

(3.3) d t u(;x)(t) = Au(t,-)(x), Nu{t,-)\ dn = 0,u{0,-) = f. 
Then 

P t f(x) = Ef(X* t ), t>o,xen. 
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Proof. Since u(t,x) := Ptf(x) satisfies (13.31) and Xf satisfies (13.11) . by the Ito formula, 
the process 



<j>(s):=u(t-s,Xl) = PJ(x)+ [ S 2(Vu(t-s,-)(Xl),dW s ), s e [0,t] 

Jo 

is a martingale. In particular, 

P t f{x) = E$(0) = E$(t) = Ef{X£). 



□ 



Proof of Theorem I i.£ As explained in the proof of Theorem II. 1\ it suffices to prove (II. ip 
for t e (0, 1]. To this end, we shall make use of the reflecting Brownian motion introduced 
above. 

Let Ptf(x) = Kf(X^). To estimate the gradient of P t f, let y ^ x be two points in 
Q, and let (X* , l*) and (X$ , If) solve (13. ip with X = x and y respectively. Then, as 
explained in the proof of Lemma 13. 1| the convexity of Q implies that 

\X?-X?\ < \x-y\, t>0. 



Thus, 



\P t °f(y) - P?m\ < E|/(xf ) - f(x?)\ < E (\f(x?) - f(xf 



\x -y\ \x- y\ V \Xf - X, 



Letting y — > x we arrive at 

(3.4) |VP°/|<i?|V/|, t>0. 
Due to an argument of Bakry- Emery [I] , this implies that 

(3.5) t\ vpo/l 2 < P t °f 2 - (if/) 2 , t>oje a b (n). 

Indeed, for a smooth function g on Q satisfying the Neumann boundary condition, (13. ip 
and the Ito formula for g(X^) imply that 

(3.6) P " g = \ f tp * A 9ds. 



o 



Since Pt = P®t, by Lemma [3.21 we have 

d 



lt P« t g = AP° t g. 
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Combining this with ( 13. 6 j) we obtain 



(3-7) ^ = \&J»g = l^Ag. 

Hence, it follows from (13. 4p and the Jensen inequality that 



(3.8) 



t|VP f °/r= [ \VP°(PU)\ 2 ds< f P a s \VPU\ 2 ds. 
Jo Jo 



(I3T71) also implies that f- ( P°{P t °_ J) 2 ) = P°| VP ; °„ jf. This, together with fETS]) , shows 



that 



Ptf ~ {Ptff = f ^P°s{Plsf) 2 ds= f "p°\VPU\ 2 ds>t\VPr 12 
Jo ds Jo 



f\ 2 - 



'0 JO 

So, (13.51) holds. Combining ( 13. 5 1) with the known uniform upper bound of the Neumann 
heat kernel on convex domains (see e.g. [HJ Theorem 3.2.9]) 

(3.9) np/iu < cr d ' 2 [ \f\(x)dx, t e (o,i],/ g & b {n), 

Jn 

we conclude that 



livp/iu = \\vp t °(p t/2 f)U<^p?[(Pt/2f) 2 ] 

< ^= t \\P t/ 2f Woo < cr^ 2 J \f\(x)dx, fea b (n) 

holds for some C > and all t G (0,1]. This implies the desired gradient estimate for 
t 6(0,1]. □ 



4 Applications 

Throughout this section, we let Q be a compact convex domain in M. d , and let Q° be the 
interior of Q. It is well known that the generator (A, ^(A)) of P t in L 2 (Q) is a negatively 
definite self-adjoint operator with discrete spectrum. Let A N = —A, which is thus a 
positive definite self-adjoint operator such that P t = e~ Ajvi , t > 0. 
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4.1 The Hardy spaces on compact convex domains 

For < p < 1 we let h p (R d ) denote the classical (local) Hardy space in R d . We consider 
two versions of this space adapted to Q. Let denote the space of C°° functions with 

compact support in Q°, and let denote its dual, the space of distributions on Q°. 

The first adaptation of the local Hardy space to Q, denoted by h p (Q), consists of elements 
of which are the restrictions to Q° of elements of h p (R d ). That is, for < p < 1 we 

set 

ftP(Q) : = {/ e y'{R d ) : there exists F G h p (R d ) such that F\ n o = /} 
= h p {R d )/{F G h p (R d ) : F = in 

which is equipped with the quasi-norm 

H/ll^o) := M{\\F\\^ m : F G h p (R d ) such that F\ Q o = /}. 

The second adaptation of the local Hardy space to Q, denoted by h p (Q), consists of 
distributions in Q° with extension by zero to R d belonging to h p (R d ). More specifically, 
for < p < 1, 

h p (Q) := h p {R d ) n {/ G h p (R d ) : / = in Q c }/{f G h p (R d ) : / = in 

where Q c := R d \ Q. We can identify h p (Q) with a set of distributions in which, 
when equipped with the natural quotient norm, becomes a subspace of h p (Q) (see [Hj for 
more details). Obviously, we have that h p (Q) C h p (Q) whenever d/(d+l) < p < 1. 

Given a function / G L 2 (f2), consider the following local version of the non-tangential 
maximal operator associated with the heat semigroup generated by the operator A at: 

N loCjAN f(x) := sup |e-* 2A -/(y)|, x e Q. 

y€Cl, \y—x\<t<l 

For < p < 1, the space ^Ajv(^) ^ s then defined as the completion of L 2 (Q) in the 
quasi-norm 

(4-1) ll/lk„(n) == l|iVioc,A N /|U, ( c). 

Based on the gradient estimate f ll.2p in Theorem 11.21 and the adapted atomic theory 
of Hardy spaces Ii p An (Q), the following result is obtained in [T2], see also [21 El EEH EE2] and 
references within for discussions on the Hardy spaces h p (Q). 

Theorem 4.1 (|12j). Let Q be a compact convex domain in R d . Then h p (Q) = h p A (Q) 
ford/(d+l) <p<l. 



4.2 Riesz transforms on bounded convex domains 

Consider the general 
cian A at, defined by 



— 1/2 

Consider the generalised Riesz transform T = VA^ associated to the Neumann Lapla- 



9 r°° /7q 
o V s 



The operator VA^ 2 is bounded on L 2 (Q). In [9], it is shown that the operator VA^ 2 
is of weak (1,1) by making use of (13.91) . hence by interpolation, is bounded on L P (Q) for 
1 < p < 2. For the case p > 2, according to [2], the following assertions are equivalent: 

(1) For all p G (2, oo), there exists C p such that 



|Ve-* Ajv ||| < 



— 1/2 

(2) The Riesz transform VA^y is bounded on L p (f2) for p e (2, oo). 



In terms of the gradient estimate (11. 2p in Theorem 11.21 it deduces the following theo- 



rem. 



Theorem 4.2. Let Vt be a compact convex domain in W 1 . Let T = VA^ be the Riesz 

— 1/2 

transform associated to the Neumann Laplacian Ajy on Q. Then the operator VA^ is 
bounded on L P (Q) for all 1 < p < oo. 

— 1/2 

If p = 1, then the operator VA N is also of weak type (1, 1). 

By (ll.2p and (13. 9p . we can also extend the Riesz transform VA^ 2 to a bounded 
operator from h p AN (Q) into L P (Q) for < p < 1. Hence by Theorem 14.11 it can be 
extended to a bounded operator from h p (Q) into L P (Q) for d/(d + 1) < p < 1. The proof 
is similar to that of [T2~| Theorem 4.2]. 



4.3 Regularity of solutions to the Neumann problem 

Define the Neumann Green operator as the solution operator C°°(Q) 3 f h- > u — 
Gat(/) £ W 1,2 {Vl) for the Neumann problem 




An = / in Vt° 
Nu =0 on dtt, 



where it is also assumed that J n f = Q and the solution is normalized by requiring that 
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4.3.1 Estimate for the gradient of Green potential 



Let W S ' P (Q) stand for the Sobolev space of functions in L P (Q) with distributional deriva- 
tives of order s in L P (Q). By L P ± (Q) and W^ P (Q) we denote the subspaces of functions / 
in L p (n) and W s ' p (tt) subject to J n f(x)dx = 0. 
Recently, V. Maz'ya proved the following result. 

Theorem 4.3 Ql5j). Let Q be a convex domain in IR rf . Let f £ L q ± (Q) with a certain 
q > d. Then there exists a constant c depending only on d and q such that the solution 
u £ W]^ 2 {yt) of the problem satisfies the estimate 



(4.3) ||V«|U» (n) < c{d,q)C^ l \^- d ^ d 



1,8(0)- 



Note that by Theorem ll.2[ we can give a simple proof of Theorem 14.31 Indeed, based 
on the gradient estimate (11. 2p of p t in Theorem 11.21 we have 

C 

\VG N (x,y)\ < _ Vx,y£fi, 
\x y\ 

where GV is the Green function for the Neumann semigroup on Q, and hence by a standard 
argument, estimate (14. 3[) follows readily. 



4.3.2 Estimate for the second-order derivatives for Green potential 

Let Gn be the Neumann Green operator for the Neumann problem in (14.21) . Note that 
L 2 -boundedness of the mappings 

r , t d 2 G N (f) 

has been known since the mid 1970's ((IHj), but optimal L p estimates, valid in the range 
1 < p < 2, have only been proved in the 1990's by Adolfsson and D. Jerison pp. It 
should be mentioned that the aforementioned LP continuity of two derivatives on Green 
potentials may fail in the class of Lipschitz domains for any p £ (1, oo) and in the class of 
convex domains for any p £ (2, oo) (see [1] for counterexamples; recall that every convex 
domain is Lipschitz). 

A natural question is to study the regularity of the Neumann Green operator when 
the L p -scale is replaced by the scale of Hardy spaces, H p , for < p < 1. Recently, X. 
Duong, S. Hofmann, D. Mitrea, M. Mitrea and the second named author of this article 
gave a solution to the conjecture made by D.-C. Chang, S. Krantz and E.M. Stein([8]) 
regarding the regularity of Green operators for the Neumann problems on h p (Q) and 
h p (Q), respectively, for all < p < 1, and this range of p's is sharp (see [T2]). 
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Theorem 4.4 (|12j). Let Q be a compact, simply connected, convex domain in M d and 
recall that G^r stands for the Green operator associated with the inhomogeneous Neumann 
problem Then the operators 



(4.4) £%L, i,j = l,...,d, 

originally defined on {/ G C°°(Q) : j n fdx = 0}, extend as bounded linear mappings 
from h"P(Q) to h^(Q) whenever d/(d + 1) <p < 1. 

/fp = 1, then the operators g x . d ^, , hj = 1,. . .,d, are also of weak type (1, 1). Hence 
by interpolation, they can be extended to bounded operators on L P (Q) for 1 < p < 2. 

The proof of this theorem is obtained by using suitable estimates for singular integrals 
with non-smooth kernels and an optimal on-diagonal heat kernel estimate. We should 
mention that the gradient estimate ( jl .2j) played a major role in the proof of this theorem, 
regarding the regularity of Green operators for the Neumann problems on hf(Jl) and 
h p z {Q), respectively, for all < p < 1. For the detail of the proof, we refer the reader to 



Acknowledgments. L.X. Yan would like to thank XT. Duong, S. Hofmann, D. Mitrea 
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